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ABSTRACT 



| I | , Context. Several young supernova remnants (SNRs) have recently been detected in the high-energy (HE; 0.1 < E < 100 GeV) and 

very-high-energy (VHE; E > 100 GeV) gamma-ray domains. As exemplified by RX J 17 1 3. 7-3946, the nature of this emission has 
been hotly debated, and direct evidence for the efficient acceleration of cosmic-ray protons at the SNR shocks still remains elusive. 
Aims. We study the broadband gamma-ray emission from one of these young SNRs, namely RCW 86, for which several observational 
^ ■ lines of evidence indirectly point towards the presence of accelerated hadrons. We then attempt to detect any putative hadronic signal 

j.^ ■ from this SNR in the available gamma-ray data, in order to assess the level of acceleration efficiency. 

Methods. We analyzed more than 40 months of data acquired by the Large Area Telescope (LAT) on-board the Fermi Gamma-Ray 
Space Telescope in the HE domain, and gathered all of the relevant multi-wavelength (from radio to VHE gamma-rays) information 
about the broadband nonthermal emission from RCW 86. For this purpose, we re-analyzed the archival X-ray data from the ASCA/Gas 
Imaging Spectrometer (GIS), the XMM-Newton/EPlC-MOS, and the ^re/Proportional Counter Array (PCA). 
Results. Beyond the expected Galactic diffuse background, no significant gamma-ray emission in the direction of RCW 86 is detected 
^ | in any of the 0.1-1, 1-10 and 10-100 GeV Fermi-LAT maps. The derived HE upper limits, together with the H.E.S.S. measurements 

l/") , in the VHE domain, are incompatible with a standard E~ 2 hadronic emission arising from proton-proton interactions, and can only 

be accommodated by a spectral index T < 1.8, i.e. a value in-between the standard (test-particle) index and the asymptotic limit 
C^l ' of theoretical particle spectra in the case of strongly modified shocks. In such a hadronic scenario, the total energy in accelerated 

\ particles is at the level of r] CR = E C r/E sn ~ 0.07 d^ 5kpc /" C m-3 (with the distance d 2 5k pc = d/2.5 kpc and the effective density fi cm _ 3 = 

■ ;1/1 cm -3 ), and the average magnetic field must be stronger than 50 pG in order to significantly suppress any leptonic contribution. 

On the other hand, the interpretation of the gamma-ray emission by inverse Compton scattering of high energy electrons reproduces 
| the multi- wavelength data using a reasonable value for the average magnetic field of 15-25 pG. In this leptonic scenario, we derive 

a conservative upper limit to t]cr of 0.04 d 2 5k /n cm -3- We discuss these results in the light of existing estimates of the magnetic field 
strength, the effective density and the acceleration efficiency in RCW 86. 



Key words, gamma rays: observations - supernova remnants - ISM: individual (RCW 86, G315.4-2.3, MSH 14-65) 

1. Introduction linear version NLDSA, see iMalkov & Drurvl l2001h . which is 

commonly invoked to explain several observational (though, in- 

Supernova remnants (SNRs) are thought to be the primary direct) lines of evidence for efficient particle acceleration at the 

sources of the bulk of Galactic cosmic-ray (CR) protons ob- SNR forward shocks up tQ yery high energies Amongthereq- 

served at Earth, up to the knee energy at ~ 3 PeV. This paradigm uiskes that must be MMed by Ms theory (fjee Bl as jggf^ for 

mainly relies on the need to have sufficiently energetic sources a recent ^ the observe d broadband nonthermal emission 

that could provide the necessary power to maintain the Galactic of individual SNRs is of particu i ar inter est. This emission, aris- 

CR energy density (e.g. | Fields et al. || 2001|) , and on the domi- ing from accelerated particles that emit photons through sev . 

nance of SNRs among sources of non-thermal radio emission. eral channels (sync hrotron [SC], inverse-Compton [IC], non- 

Our understanding of CR acceleration in SNRs mainly relies on thermal bremsstrahlung> pro ton-proton interactions and subse- 

the so-called Diffusive Shock Acceleration theory (in its non- quent ^ decayX offers new insights into the partide acceleration 

mechanisms at work in these sources, given the large number 

Send offprint requests to: lemoine@cenbg.in2p3.fr, 
mrenaud@lupm.univ-montp2.fr 
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of multi-wavelength observations available nowadays towards 
many Galactic SNRs. 

In particular, recent observations of young SNRs in the 
high-energy (HE; 0.1 < E < 100 GeV) and very-high-energy 
(VHE; E > 100 GeV) gamma-ray domains have raised sev- 
eral q uestions and triggered numerous theoretical investigations 



eral questions and triggered numerous theoretical investigations 
(e.g. IZirakashvili & Aharonianll2010t iBerezhko & VoUdl2010t 
iFang et al.ll2.01 lb ICapriolill20 1 lh . The critical issue regarding the 
nature of the observed gamma-ray emission has been intensively 
discussed in the literature in recent years (Ellison e t alJ l2.Q10l 
and references therein). Nevertheless, these joint HE/VHE ob- 
servations can help us to discriminate spectrally between the 
leptonic (through IC emission) and had ronic (through n° — >2y) 
scenarios, as shown by the Ferm z-LAT (lAbdo et al. 1201 lal) and 
H.E.S.S. ( lAharonian et alj|2007l) observations towards the TeV- 
bright SNR RX J17 13.7-3946, which tend t o support a leptonic 
model for the observed gamma-ray emission dEllison et al.l2012t 

iLi et alj[20TTh . 

RCW 86 dRodgers et al.l 19601) . also known as G3 15.4-2.3 or 
MSH 14-65 (iMills et al.ll 19611) . is a 42' diameter Galactic SNR 
in the southern sky. There has been much controversy about the 
nature of the supernova (SN) progenitor, and the SNR age and 
distance. This has resulted from the difficulties in reconciling the 
young age of RCW 86, based on its connection with the first his- 
torical SN ever recorded in AD 185, with its large size, given 
the relatively large distance estimates, which place RC W 86 at 
2.3-2.8 kpc dRosado et al.ll 1996t ISollerman et al.ll2003l) near an 
OB stellar association. IWilliams et al.l (1201 ll) recently reviewed 
all the arguments about the nature of the SN progenitor and criti- 
cally examined the available observations of RCW 86. They sug- 
gest that the SNR is likely the remnant of a type la SN, and, 
from hydrodynamic simulations, that the off-center explosion 
occurred in a low-density c avity carved by the progenitor system 
(see also IVink et al.lfl997l) . This scenario allows one to explain 
at once the young age and the large distance of RCW 86, which 
we scale in terms of d2.5 = d/(2.5 kpc) throughout the paper. 

The g eneral outline of its nearly circular shell is similar in 
the radio (IWhiteoak & Green 199 6t lDickel et al.ll2001l) . infrared 
(IR, IWilliams et al.ll201 lh, optical (e.g. 



rav rtVink et alJI 1997b iBamba et al 



Borkows ki et al.ll2001bb iRho et al 



2000; 



2002 



1 et al.ll200lh . 
Srru^lT997l) 



and X- 



Bocchino et al]l2000t 



Vink et al.! 2006) do 



mains, alth ough significant fine-scale differences have been re- 
ported (e.g. IRho et al.ll2002l and Fig.©. X -ray observations to- 
wards RCW 86 have revealed the presence of both thermal and 
nonthermal emission, with very distinct morphologies: while 
soft X-rays correlate with optical emission from non-radiative 
shocks and IR emission from collisionally heated dust, the con- 
tinuum hard X-ray emission is seen at the edges of radio emis- 
sion. The high-temperature plasma, which mostly contains the 
strong Fe K a line emission, also shows a particular morphology 
extending towards the SNR interior, as revealed with Suzaku, and 
is thou ght to originate from Fe-rich ejecta heate d by the reverse 
shock dUeno et al .120071; I Yamaguchi et al.l2008h . The global dis- 
tribution of the Fe-rich pl asma in the entire SNR has recently 
been mapped with Suzaku (I Yamaguchi et alj201 ll) . and the total 
Fe mass of ~ 1 M Q deduced from these observations, together 
with a relatively low ejecta mass of 1-2 M Q , strengthens the sce- 
nario of a type la SN at the origin of RCW 86. 

Physical conditions (shock speed, ambient density, and mag- 
netic field) vary greatly along the shell-lik e structure. In partic- 
ular, slow shocks (~ 6 00-800 km s -1 , see lLong & Blairll 1990b 
avamian et al. 2 001b and relative ly high post-shock densities 



|Gh 



Williams et al. I I20TT1) ' have been measured in the 



east (NE) region exhibits much faster shocks (> 2700 km s 1 and 
6000 ± 2800 km s _1 , see IVink et all 120061: IHelder et al.ll2009l) 
and l ower densities ( ~ 0.1-0.3 cm 3 , e. g. lYamaguchi et al.l 
120081) . In this region, IHelder et all d2009l) argued that at least 
50% of the total pressure is induced by CRs, based on the dis- 
crepancy between the measured shock velocity and the spectro- 
scopically determ ined post-s hock proto n temp erature (see also 
IViriket all 1201 Oh . Moreover. IVink et al.l (120061) found evidence 
for a concave electron spectrum, as predicted by the NLDSA 
theory, in order to explain the radio and X-ray SC emission ob- 
served from the same region. These two measurements, together 
with the recent detection of RCW 86 with the H .E.S.S. experi- 
ment in the VHE domain (lAharonian et al 1 l2009h . seem to point 
towards an efficient CR source. However, complementary HE 
observations are needed to probe the nature of the gamma-ray 
emission, as discussed above. 

We here report on Fermz'-LAT observations and data analy- 
sis towards RCW 86 (section In an attempt to constrain the 
broadband nonthermal emission from the SNR, we present a re- 
analysis of the archival ASCA/GIS, XMM-Newton/EPlC-MOS, 
and RXTE/PCA X-ray data (section [3}, and collect the available 
information in the radio and VHE gamma-ray domains. We then 
discuss the constraints derived on the parent particle spectrum 
and on the subsequent acceleration efficiency, in the light of ex- 
isting estimates (section[4]i. 



2. Fermi-LPJ observations and data analysis 

2.1. Observations 

The LAT is a gamma-ray telescope that detects photons by con- 
version into electron-positron pairs in the energ y range from 
20 Me V to more than 300 GeV, as described by lAtwood et al] 
(I2009I) . It is made of a high-resolution converter/tracker (for 
direction measurement of the incident gamma-rays), a CsI(Tl) 
crystal calorimeter (for energy measurement) and an anti- 
coincidence detector to reject the background of charged par- 
ticles. The LAT has a large effective area (~ 8000 cm 2 on-axis 
above 1 GeV), a wide field of view (~ 2.4 sr) and good angu- 
lar resolution (~0.6° radius for 68% containment at 1 GeV for 
events converting in the front section of the track er). The on- 
orbit calibration is described in LAbdo etalJd2009ah . 

The analysis used more than 40 months of data collected 
starting on August 4, 2008, and extending until February 22, 
2012. We selected events with energies greater than 0. 1 GeV, and 
excluded those with zenith angles larger than 100°, to minimize 
contaminatio n from secondary gamma rays from the Earth's 
atmosphere (lAbdo et al 1 l2009bl) . In addition, we excluded the 
time intervals when the rocking angle was more than 52° and 
when the Fermi satellite was within the South Atlantic Anomaly. 
We used the P7SOURCEJV6 Instrument Response Functions 
(IRFs), and selected the 'Source' events, which correspond to 
the best compromise between the number of selected photons 
and the charged particle residual background for the study of 
point-like or slightly extended sources. 



2.2. Analysis and results 

For the spectral analysis, a spectral-spatial model containing dif- 
fuse and point-like sources was created, and the parameters were 
obtained from a maximum likelihood fit to the data using the gt- 
like code implemented in the Fermi Science Toolfl The like- 



(~ 2 cm , see 

southwest (SW) and northwest (NW) regions, while the north- 1 http://fermi.gsfc.nasa.gov/ssc/data/ 
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lihood tool gtlike fits a source model to the data along with 
models for the residual charged particles and diffuse gamma-ray 
emission. For the Galactic diffuse emission, we used the ring- 
hybrid model galJ2yearp7\>6_v0.fits. The instrumental back- 
ground and the extragalactic radiation are described by a sin- 
gle isotropic component with the spectral shape in the tabu- 
lated model iso-p7v6source.txt. These models are available from 
the Fermi Science Support Centefl All sources within 15° of 
RCW 86 were ex tracted from the Fermz'-LAT 2FGL catalog 
dNolan et al.ll2012l) and added to our spectral-spatial model of 
the region. The spectral parameters of sources closer than 3° to 
RCW 86 were allowed to vary in the likelihood fit, while the 
parameters of all other sources were fixed at the 2FGL values. 
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Fig. 1. Fermi-LAT test statistic (TS) map towards RCW 86, 
above 1 GeV. The TS was evaluated by placing a point-source 
at the center of each pixel, Galactic diffuse emission and 
nearby 2FGL sour ces being included in th e background model. 
H.E.S.S. contours (lAharonian et al.ll2009l) are overlaid as white 
solid lines, and the white circle indicates the size of the 68% 
containment region of the LAT PSF at 10 GeV. Crosses and el- 
lipses mark the best-fit positions and associat ed errors of the two 
sources from the 2FGL Fermi-LAT catalog dNolan et al.ll2012|) 
within 4° from RCW 86. 



Using this spectral-spatial model, we computed the LAT test 
statistic (TS) map of the region containing RCW 86, as shown 
in Figure Q] The TS is defined as twice the difference between 
the log-likelihood L\ obtained by fitting a source model plus 
the background model to the data, and the log-likelihood Lq ob- 
tained by fitting the background model only, i.e. TS = 2(Li - 
Lq). A TS = 25 with two degrees of freedom (spectral index and 
normalization) roughly corresponds to a 4.6cx detection signif- 
icance. Figure [T] contains the TS value for a point source with 
fixed photon index F = 2 at each map location, thus giving 
a measure of the statistical significance for the detection of a 
gamma-ray source in excess of the background. No significant 
excess suggesting gamma-ray emission is evident in the skymap, 
apart from faint excesses, which are unassociated with any 2FGL 
catalogued source, at the northern and southern borders of the 
map, and a faint hotspot located near the NE shell of RCW 86. 
Under the assumption of a point-like source, this excess is found 
at R.A.(J2000) « 221.0° and Dec(J2000) ~ -62.4° with a TS 



value of 12 {i.e. ~ 2.4 <x for four dof, where the spectral index 
is allowed to vary). Therefore, we do not consider this hotspot 
in the following, bearing in mind that further Fermi-LAT obser- 
vations in the coming years might provide confirmation of this 
indication of a signal towards RCW 86. 

Since the Fermz'-LAT 2FGL catalog does not contain any 
source in a region of 1 around the SNR RCW 86, to determine 
an upper limit to the gamma-ray emission from RCW 86, we 
added an extended source consistent with the H.E.S.S. observa- 
tions, i.e. a uniform disk centered at (R.A.(J2000) = 220.68°, 
Dec.(J2000) = -62.44°) with a radius of 0.40°. We modeled the 
emission from RCW 86 using a power-law spectral distribution 
with two free parameters: the flux and the spectral index. The 
test statistic was employed to evaluate the significance of the 
gamma-ray fluxes coming from RCW 86 in three different en- 
ergy bands: 0.1 GeV - 1 GeV, 1 GeV - 10 GeV, and 10 GeV 
- 100 GeV. In all tested energy bands, the likelihood ratio test 
indicate values of TS less than 25. This agrees with the non- 
detection in the 2FGL catalog as well as the LAT TS map pre- 
sented in Figure Q] Assuming a power-law shape for the source 
with a fixed spectral index of 2, we then derived 99.9% con- 
fidence level (CL) upper limits to the flux in the se three energ y 
intervals using the Bayesian method proposed bv lHelenel J 1983). 
These upper limits, shown in Figure|4] amount to 1 .3 x 10~ 8 , 1 . 1 



x 10 , and 3.4 x 10 ph cm s , respectively. Varying the 
spectral index between 1.5 and 3 did not change significantly 
these estimates. 



3. X-ray observations 

We seek to estimate the total nonthermal X-ray emission 
from RCW 86 in order to constrain the underlying acceler- 
ation mechanisms, together with the gamma-ray observations 
presented above. There is a large amount of X-ray data to- 
wards the SNR, since ASCA, XMM-Newton, and more re- 
cently Suzaku hav e performed a nearly complete coverag e 
(iBamba et al.l 120001: IVink et al.l 120061 lYamaguchi et al.l 1201 ll) . 
In hard X-rays, the non-imaging instrument PCA onboard 
RXTE detected R CW 86 abov e 10 ke V, and the original results 
were presented in lAllen et al ] (119991) . From our point of view, 
it is worthwhile revisiting these archival X-ray data for sev- 
eral reasons. The RXTE/PCA data analysis presented below was 
carried out with the latest software release, and we investigate 
(and hence correct the extracted spectrum for) the contamina- 
tion from the Galactic Ridge and the effect of the PCA off-axis 
response for such an extended source. For both ASCA/GIS and 
RXTE/PCA datasets, we also want to constrain the SC emission 
with the most advanced, phy sically-motivated, modeling dedi - 
cated to X-ray observations (Zir akashvili & Aharonianl 1201 Oh . 
rather than either the usual power-law approximation, which 
is obviously i nadequate in the cutoff regim e, or the often-used 
srcut model dRevnolds & Keohandfl999l) . which relies on the 
SNR properties in the radio domain and hence makes the im- 
plicit assumption that radio and nonthermal X-rays arise from 
the same emission region , at odds with what is observed in 
RCW 86 dRho et al.ll2002l and Fig.©. 



3.1. ASCA/GIS 



RCW 86 was observed with ASCA dTanaka et alj|1994l) on 1993 
August 17-18 with three pointings towards the NE, western (W), 
and SW regions. ASCA has two kinds of in struments, GIS (Gas 
Imaging Spectrometers: IOhashi et al] 19961) and SIS (Solid-state 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModtffltaguag Spectrometers: Burke et al.ll 199 lh . each at the focus of 
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an X-ray Telescope dSerlemitsos et al.lll995l) . We focus on the 
GIS data set since it has a larger field of view and effective area. 
Data selection and analysis were perfo rmed following the pro- 
cedure described in lBamba et al.l feOOO). The resulting exposure 
times are -20 ks for NE, ~9 ks for W, and ~8 ks for S W. We pro- 
duced ASCA/GIS spectra from three distinct regions of RCW 86, 
which we labeled SW, NE and inner, as shown in Figure [2] 
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Fig. 2 . Composite radio (MGPS-2 at 843 MHz, iMurphv et ail 
I20071 in red) and X-ray (XMM-Newton/MOS in the 2.5-4.5 
keV band, in blue) image of RCW 86. The XMM-Newton/MOS 
image has been smoothed with a Gau ssian of width 10". 
H.E.S.S. contours (lAharonian et al 1 120091) are overlaid in white, 
as in FigureQ] Triangles mark the RXTE pointing positions with 
the associated exposure times. The different ASCA/GIS regions 
for spectral extraction are delineated in green and labeled NE, 
SW and Inner (counts w ithin the dashed circle were excluded), 
as in lBamba et ai1(l2000h . 



Since RCW 86 exhibits both thermal and nonthermal X- 
ray emission at different relative levels along the shell, we fit 
these spectra individually with an absorbed three-component 
model consisting of (1) a plane -parallel shock model with vari- 
able metal abundanc es (with their relative a bundances fixed to 
solar ratios, pshock. fBorkowski et al.ll2001al) . (2) an unresolved 
Gaussian line at fixed energy of 6.4 keV, and (3) the analytical 
appro ximation of the SC emission (IZirakashvili & Aharonianl 
1201 Ol ZA10 her eafter) from the shock-acce le rated electron spec- 
tra derived by IZirakashvili & Aharonianl d2007l) in the SC- 
loss-limited regime. The first two components account for the 
low-temperature thermal emission in non-equilibrium ioniza- 
tion and th e Fe-K line detected in several regions a long the 
SNR shell dBorkowski et al.l l2001bt lUeno et al.lT2007l) . respec- 
tively. We note that this iron line is m ost likely produced 
by a high-temperature (> 5 keV) plasma (Bamba et al1l2000t 
iBocchino etlri1l2000t IBorkowski et alJl2001bl iRho et aljl2002l) . 
which we do not model here. However, we checked that adding 
such a component does not change significantly the best-fit re- 
sults. The third component was used t o describe the X-ray SC 
emission measur ed in RCW 86 (e.g. IBorkowski et alj 12001 tit 
IVink et al.ll2006l) . We followed the ZA10 prescription, accord- 



ing to which the SC spectrum from shock-accelerated electron s 
in the downstream region (see Zirakas hvili & Aharonianll2007l) . 
in the case of Bohm diffusion and for a distributed (turbulent) 
magnetic field, is found to have the following behavior: 



dN(E) 
dE 





\ (Ef 4 ~ 


25/4.8 




oc E~ 2 


1 + 0.185(-) 


x exp 





Ei = 



0.056 keV 



>1B 



V 



f 



3 x 10 3 km s~ 



(1) 
(2) 



for which the upstream to downstream magnetic field ratio k is 
assumed to be 1/ VTT, Vf is the forward shock velocity, and tjb = 
D/Ds hm > 1 is the factor allowing for deviation of the diffusion 
coefficient from its nominal value. The above expression is valid 
when the electron spectrum in the cutoff region is determined 
by the energy losses, i.e. when the SNR age is older than the 
SC cooling time of electrons t ~ 1300 rj l J 2 (V//10 3 km s^ 1 )" 1 
(B/25 fiGy 3 ^ 2 yr. Provided that the diffusion proceeds close to 
the nominal Bohm regime, this condition is met for RCW 86, 
which exhibits shock speed s in excess of 700 km s -1 and B-fiel d 
estimates of at least 25 fiG (IVolk et. al.ll2005tlVink et~aT]l2006h . 

The ASCA/GIS spectra and their respective best-fit models 
are shown in Figure[3] while the spectral parameters are provided 
in Table [T] Results are in good agreement with tho se reported 
bv lBambaetall (l2000l) and Borkowski et all (l2001bl) . who both 
used ASCA data and different spectral models. Regarding the 
nonthermal emission, the best-fit values of Ei ~ 0.01-0.02 keV 
translate into shock speeds on the order of 1300-1800 -s/tJb km 
s _1 , which are roughly compatible with previous estimates for 
the different parts of the RCW 86 shell. As shown in Table Q] 
the ionization ages r from the pshock models are poorly con- 
strained, and we fixed them to their nominal values while ex- 
tracting the errors in the other free parameters. Finally, the three 
defined regions do actually not encompass the whole SNR. To 
estimate the nonthermal flux outside the ASCA regions, we used 
the 2.5-4.5 keV XMM-Newton/MOS image shown in Figure [2] 
assuming that the flux in this band is entirely nonthermal. We 
found that 80% of the 2.5^1.5 keV flux is contained within the 
ASCA regions, and we consider in the following the best-fit un- 
folded spectra shown in Figure [3] (bottom) corrected by this fac- 
tor. 

3.2. XMM-Newton/MOS 

XMM-Newton observed RCW 86 during several cycles, as part 
of the guaranteed time (GT, PI: Bleeker) program in 2000 and 
2001 and as part of the guest observer program (PI: Vink) in 
2004 and 2007. The advantage of XMM-Newton observations 
over the ASCA observations is that almost the entire SNR was 
covered, and that XMM-Newton has a larger effective area. This 
allowed us to estimate the total nonthermal X-ray flux, includ- 
ing that from regions of faint nonthermal emission outside the 
ASCA regions. Unfortunately, the very elliptical and high orbit 
of XMM-Newton results in a higher background rate than ASCA, 
which, moreover, is very irregular. To obtain a nonthermal flux 
estimate as accurately as possible we limited ourselves to the 
observation IDs the least plagued by a high background sig- 
nal. The data reduction was carried out with the standard XMM- 
Newton software SAS version 20110223. To limit the effects of 
chip-gaps and variable background, we limited the analysis to 
the EPIC-MOS instruments. Since the XMM-Newton observa- 
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tions were spread out over a large time interval, with many dif- 
ferent, overlapping pointings, we based our nonthermal flux esti- 
mate on imaging rather than on a full spectral analysis, since for 
the latter one would need to take into account the different back- 
grounds per spectrum and different exposure weights for each 
single extraction region. The image for the EPIC-MOS 1 and 2 
instruments, shown in Figure [2] was made in the 2.5^.5 keV 
energy band. This band was used to estimate the total nonther- 
mal flux, which starts dominating above ~ 2 keV although some 
S XV line emission can still be present at 2.45 keV. The cut-off 
at 4.5 keV was used for practical reasons, as the relative back- 
ground signal increases at higher energies. The effective area be- 
tween 2.5 and 4.5 keV is relatively flat and we adopt here the 
mean value of 345 cm 2 . The background estimates were based 
on the regions of the image outside RCW 86. This is not ideal, 
as the background rate varies from observation to observation, 
with an estimated background level of (0.75-1.5) x 1CT 4 events 
s _1 arcsec~ 2 . The total count rate in the 2.5-4.5 keV band for the 
region containing RCW 86 is 3.4 counts s (for a total number 
of events of 3.7 x 10 5 ), which reduces to 1.2-2.3 counts s _1 after 
background subtraction. To convert this count rate to nonthermal 
flux, we used two approaches. One was based on the analysis o f 
the nonthermal emission from the NE region (IVink et al.ll2006l) . 
For the MOS spectra of this region, we have both a flux esti- 
mate and a total 2.5-4.5 keV count rate. Using the best-fit spec- 
tral index of T = 2.8, we found that the conversion factor from 
count rate to nonthermal flux in the 2.5-^.5 keV band is ~ 1.7 x 
10~ n erg cirT 2 s~7( cn t/ s )- Hence, the 2.5^4.5 keV nonthermal 
flux amounts to (2.1-3.9) x 10~ n erg cirT 2 s , corresponding 
to a normalization at 1 keV of (5.7-11) x 10~ 2 cirT 2 s keV -1 . 
One can also directly convert the net count rate to a flux rate 
using the effective area given above, which gives a normaliza- 
tion at 1 keV of (5.1-9.8) x 1(T 2 cirT 2 s" 1 keV" 1 . These values 
are consistent with the measurements by ASCA (and RXTE, see 
below and Table Q]). 
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Fig. 3. 0.5-10 keV ASCA/GIS spectra from the three regions 
considered in this work, together with the RXTE/PCA spectrum 
in the 10-25 keV energy range. All spectra have been rebinned to 
10 cr per bin for sake of clarity. Solid lines represent the best-fit 
models (see text and Table [TJ. Residuals and unfolded nonther- 
mal spectra are shown in the middle and lower panels, respec- 
tively. 



3.3. RXTE/PCA 

The SNR was observed with RXTE for a total of 4 2 ks between 
1997 February 26 and 1997 March 17. The PCA dJahoda et al.l 
120061) is comprised of an array of five coaligned proportional 
counter units (PCUs) that are mechanically collimated to have 
a field of view of 1° (FWHM). For RCW 86, the on-axis effec- 
tive area has a peak of about 4600 cm 2 at 7.3 keV and is greater 
than or equal to 10% of this value from about 2.4 to 27 keV. 
Between these energies, the fractional energy resolution AE/E 
ranges from about 0.25 to 0.12 (FWHM), respectively. The data 
for observation ID 20260 were downloaded from the HEAS ARC 
archive and reprocessed using version 6.11 of FTOOLSE The 
tools xtefilt and raaketime were used to produce a new filter 
file for each one of the five pointings in the observation (whose 
locations are shown in Figure [2) and a set of good-time inter- 
vals (GTIs) for standard event selections, respectively. The GTIs, 
which include 38 ks of data, and the tool saextrct were used to 
obtain a single, combined spectrum for the first layer of anodes 
in all five PCUs. An estimate of the combined spectrum of the 
instrumental and cosmic X-ray backgrounds was produced us- 
ing the tool pcabackest with the latest combined faint-source 
model and SAA history file. As a test of the accuracy of this 
model, the total number of background counts in the last 50 
of 129 pulse-height channels (i.e. at energies above 43 keV) 
is 698,135. For comparison, the total number of counts in the 



3 http: //heasarc. nasa.gov/ftools/ftoolsjnenu. html 



data in the same range is 697,541. The two values agree to better 
than 0. 1 % {i.e. to better than 1 cr). A response file was produced 
using the tool pcarsp. 

The spectral analysis was performed using the energy range 
from 10 to 25 keV. Some caution should be used when inter- 
preting the results of the PCA spectral analysis because: (1) 
There are spatial variations in the cosmic X-ray background. In 
the 10-20 keV energy band, these variations are 0.2 counts s 
(1 cr) for the top layer of anodes in all five PCUfl For compar- 
ison, the background-subtracted emission from RCW 86 is 1.1 
counts s _I in this energy band. If an additional power- law com- 
ponent with a photon index of 1.3 (lAiello et al.ll2008h is added 
to the fit to represent excess cosmic X-ray background emis- 
sion, then the best fit normalization of this component is essen- 
tially zero and the 90% upper limit to the normalization corre- 
sponds to 0.3 counts s _1 . (2) The background model does not in- 
clude emission from the Galactic Ridge. Unfortunately, the anal- 
ysis of the Galactic Ridge emission, as measured using the PCA 
(IValinia & MarshalHl 1998b . does not extend to the two pointing 
locations for RCW 86 (R.A.(J2000) = 220.75°, Dec.(J2000) = 
-62.46° for the center and R.A.(J2000) = 220.25°, Dec.(J2000) 
= -62.67° for the SW, see Figure However, there is a PCA 
observation (ID 30267) of a nearby, more or less source-free, re- 
gion of the sky (R.A.(J2000) = 214.55°, Dec.(J2000) = -63.93°). 

4 http://universe.gsfc.nasa.gov/xrays/programs/rxte/pca/doc/bkg/bkg-2002 
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Table 1. Spectral models for RCW 86. Uncertainties quoted here are provided at the 90% confidence level. 



An analysis of the data for this observation reveals little evidence 
of Galactic Ridge emission in the 10-25 keV energy band, and 
we assume that we can ignore emission from the Galactic Ridge 
in our analysis. (3) While RCW 86 is entirely within the field 
of view of the telescope for each pointing, the SNR is not a 
point source and the relative PCA detection efficiency declines 
more or less linearly from o ne for an on-axis so urce to zero for 
a source that is 1° off axis dJahoda et al.ll2QQ6h . Therefore, the 
amount of emission detected is less than the amount of emis- 
sion that was produced by the SNR. On the basis of the XMM- 
Newton 2.5-4.5 keV image and the RXTE pointing directions 
shown in Figure [2] and the P CA off-axis respons e in the form 
f = 1 - r/r with r = 0.965° dJahoda et al.ll2006l) . we estimate 
this reduction factor to be ~ 0.8. The RXTE/PCA spectrum is 
shown in Figure [3] and the best-fit parameters from the ZA10 
prescription are provided in Table [TJ Fitting the spectrum with 
a power-law gives acc eptable and consis tent results with those 
originally reported bv I Allen et all d 19991) : T = 3.27 ± 0.14 and 
a normalization at 1 keV of 0.21*2q? crrT 2 s keV -1 (^ 2 /dof = 
15/31). As for the total nonthermal ASCA spectrum, we consider 
in the following the best-fit unfolded PCA spectrum shown in 
Figure [3] (bottom) corrected from this reduction factor. 

4. Discussion 

To constrain the nature of the gamma-ray emission from 
RCW 86, we have collected all the available information in the 
radio and VHE gamma-ray domains, in addition to the Fermi- 
LAT and X-ray data analyses presented in the previous sections. 
Observations towards RCW 86 have been carried out with sev- 
eral radiotelescopes. The flux densities measured with MOST at 
408 MHz (86 Jy) and Par kes at 5 GHz (18.2 Jy) are extracted 
from lCaswell et all d 1975b . with respective errors of 10% (D. A. 
Green 201 1, private communicati on). The lower limits of 22 and 
28 Jy were derived at 843 M Hz dWhiteoak & Greenl[l996l) and 
1.34 GHz dDickel et al.ll2.00lh . respectively. Spectral data points 
in the VHE domain we re obtained with the H.E.S.S. experiment 
dAharonian et alj|2009l) . The spectrum is well described with a 
power-law between 0.6 and 60 TeV, with a photon index of 2.54 
+ 0.12 and a normalization at 1 TeV of (3.72 + 0.50) x 10~ 12 
cirT 2 s _1 TeV~' (statistical errors only). 

4.1. Modeling the broadband emission of RCW 86 

Two mechanisms are commonly proposed to explain the gamma- 
ray emission in young SNRs: IC scattering off ambient pho- 



tons (called leptonic scenario) and proton-proton interactions 
(called hadronic scenario). On the basis of a phenomenologi- 
cal approach, we aim to constrain the average magnetic field 
and the total energy in accelerated particles responsible for the 
broadband nonthermal spectrum of RCW 86 shown in Figure [4] 
For this purpose, we modeled the 7r°-decay gamma-ray emission 
from inelastic interactions of accelerated hadrons with ambient 
gas nu clei, according to the procedure described in iHuang et al.l 
(I2007h . The authors calculated the gamma-ray emission from the 
collisions between CR protons and Helium nuclei and the inter- 
stellar medium with standard composition. We also computed 
the exact expressions of SC and IC emissions from a pare nt elec- 
tron spectrum, as given in iBlumenthal & Gould! (fl970l) . since 
the contribution to the IC scattering from the local interstellar 
(optical and infrared) radiation fields (ISRF), at the location of 
RCW 86, is unknown, we considered two cases: one with the 
Cosmic Microwave Background (CMB) as the only source of 
seed photons for IC scattering, and the other with both the CMB 
and the Galactic ISRFs provided by iPorter & Strond d2005l) . 
which comprise the infrared (from dust, at T ~ 35 K) and opti- 
cal (from stars, at T ~ 4600 K) emissions with respective energy 
densities of 0.66 and 0.94 eV cm~ 3 . Nonthermal bremsstrahlung 
is neglected here given the relatively low densities (< 1 cm -3 ). 

We qualitatively fit the free parameters of these leptonic and 
hadronic scenarios to the available data on RCW 86, and discuss 
the corresponding constraints on the parent particle spectrum, 
energetics, and magnetic field, in terms of the existing estimates. 
The particle spectra are assumed to follow a power-law with an 
exponential cutoff dN/dE oc E~ r x exp(-E/E max ), with the same 
spectral index for both electrons and hadrons. We note that, in 
such a case, the resulting SC emission spectrum follows a power- 
law with an exponential cutoff in the form exp(-(£'/£ , cut ) 1 '' 3 ), 
which is similar to the ZA 10 prescription that we used in section 
|3]to fit the ASCA and RXTE spectra. We denote the fraction of the 
total energy injected into accelerated particles rf ep = W etP /EsN, 
where E$n is the energy of the SN explosion, assumed to be 
the canonical value of 10 51 er|0 We also calculate the so-called 
electron-to-proton ratio K ep at momentum 1 GeV c~', which is 
compared to the value measured in cosmic -rays (K ep ~ 10~ 2 ). 
In these hadronic and leptonic scenarios, we consider two dif- 
ferent models. The first is a one-zone model that assumes that 
radio, X-rays, GeV, and TeV gamma rays arise from the same 
emission region with a constant magnetic field. Since this model 

5 This assumption would be fairly valid in t he case of RCW 86 is 
indeed the remnant of a type la SN ( Williams et alj 1201 ll : 
lYamaguchietalJl201ll) 
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cannot account for the morphological difference between radio 
and X-rays (see Figure [2]), we also consider a two-zone model, 
which exhibits two populations of radio- and X-ray -/gamma-ray- 
emitting particles. The best-fit parameters for these two pairs of 
scenarios are given in Table|2] 

4.2. Leptonic scenario 
4.2.1. One-zone model 

We assume that SC and IC photons are produced by high-energy 
electrons confined in the same emitting region. Figure [4] (top, 
blue lines) shows that the broadband emission of RCW 86 can 
be easily accomodated using an injection spectrum with spec- 
tral index T = 2.3, in agreement with the radio measurements, 
an exponential cutoff at 20-25 TeV, and an average magnetic 
field of 15-25 fiG, depending on the level of ISRFs as seed pho- 
tons in the calculation of IC emission. This one-zone model also 
provides an upper limit to the energy injected into accelerated 
hadrons of ~ 2 x 10 49 d^ 5 /« cm _3 erg, with h the effective den- 
sity as depicted by the dash-dotted blue line in Figure |4] (top). 
This estimate is only slightly sensitive to the assumed E max p , 
which we conservatively fixed here to 100 TeV. The resulting 
electron-to-proton ratio K ep amounts to > 0.05-0. 1 . 



4.2.2. Two-zone model 

The above estimates rely on the implicit assumption that radio, 
X-ray and gamma-r ay emission regio ns occupy the same vol- 
ume. As shown by iRho et all d2002l) and visible in Figure [2] 
this does not seem to be the case for RCW 86. Therefore, we 
consider two different populations of radio- and X-ray/gamma- 
ray-emitting leptons to reproduce the broadband emission of 
RCW 86. In this scenario, we use a harder lepton spectrum to 
explain the X-ray and VHE observations, with T = 2, and a sep- 
arate electron population to account for the radio emission, with 
T = 2.2. For simplicity, we do not model any putative population 
of accelerated hadrons associated with the latter leptonic popu- 
lation, given the poor constraints obtained (see below). Figure|4] 
(top, red lines) shows a reasonable fit of the X-ray and gamma- 
ray data assuming an exponential cut-off at 20-25 TeV, and an 
average magnetic field of 15-25 //G, which is similar to that of 
the one-zone model. The direct consequence of this two-zone 
model is that it slightly relaxes the constraint on the energy in- 
jected into accelerated particles, owing to a harder spectrum. In 
this case, the estimated maximal energy injected into accelerated 
hadrons is ~ 4 x 10 49 d2 5 /« C m-3 erg, i.e. twice as large as in the 
one-zone model, and the resulting K ep > ( 3-6) x 10~ 3 get s closer 
to the observed value in CRs. Recently, IWilliams et alj d201 lb 
performed a combined IR/X-ray analysis and derived post-shock 
densities in the NW and SW regions of the SNR of 2.0 cirT 3 and 
2.4 cm -3 . If applicable to the whole SNR, these values would 
imply a maximal energy i njected into hadrons o f ~ 2 x 10 49 d 2 5 
erg, much lower than what Helder et al ] (I2009h estimated in the 
NE region, where they evaluated that > 50% of the post-shock 
pressure is produced by CRs (see below, section 0. However, 
the effective density in this region is known to be much lower 



than th e western part of the shell. For instance. I Yamaguchi et al.l 
(I2008h derived a post-shock density from the (weak) thermal X- 
ray emission of 0.26 f~ cm -3 (/ is t he fill ing factor of this 
thermal component), and iHelder et alJ d2009|) ass umed a den- 
sity of 0.1 cm -3 , as derived in Vink et al.l d2006l) . As for the 
radio-SC-emitting electrons, we found that the magnetic field 
in this zone must be stronger than 10-25 jiG (depending on the 
level of G alactic ISRFs). Th is estimate agrees with the value de- 
duced by IVink et all d2006l) from the width of the filament in 
the NE r egion (24 ± 5 uG), b ut barely agrees with the one ob- 
tained bvlBamba et al.l 120051) for the SW filament (~ 4-12 fiG). 
IVolk et. all d2005l) estimated a stronger magnetic field strength 



6 For Bohm diffusion and a compression ratio of 4, k = B„ p /B f ; „„ = 

Ddown/Dup = V VTT, as defined in Equation [2] with D [town (resp. D up ) 
the diffusion coefficient in the downstream (resp. upstream) medium. 
Given that the residence time of particles t scales as 1/D, h =(n up xt up + 
ndownXt d0 wn)/(tup+tdown) ~ 0.8 n down . Throughout the discussion, we then 
assume that h equals the post-shock density. 



(99^ 4 g yi/G) from the thickness of th e same S W filamentary struc- 
ture observed bv lRho eta II (120021) . However, as noted by these 
authors, this estimate could be reduced by a factor of ~ 2, due 
to the non-spherical geometry, and hence the incomplete de- 
projection, of this featu re, which lies in the in terior of the rem- 
nant. We also note that lArbutina et al.l d2012l) found a volume- 
averaged magnetic field in RCW 86 of ~ 70 fiG, based on mod- 
ified equipartition calculations applied to shell-type SNRs. 



4.3. Hadronic scenario 
4.3.1. One-zone model 

In the hadronic scenario, gamma rays are predominantly ra- 
diated through decays of n° mesons produced in collisions 
between accelerated protons/nuclei and ambient gas. This di- 
rectly implies that the average magnetic field needs to be much 
stronger than the previous estimates to suppress the IC compo- 
nent. In a first step, as performed above, we assume that high- 
energy electrons and protons are confined in the same emit- 
ting region. The standard E~ 2 -like spectrum with a fraction 
of total energy into hadrons of r\ p ~ 0.15 d^ 5 /« cm _3, which 
would accommodate th e H.E.S.S. measurements (as discussed in 
lAharonian et alj|2009l) . is clearly ruled out, as shown in Figure 
[4] (bottom, dotted black line). The hard spectral index required 
is in-between the standard (test-particle) index and the asymp- 
totic li mit of T = 1 .5 predicted by extremely effic ient CR accel- 
eration dMalkovll 19991 iBerezhko & Ellisonlll999h . This is sim- 
ilar to RX J1713.7-3946 dAbdoetalJl2011al) an d, to a lesser 
exten t, to RX J0 8 52.0- 4622 (aka Vela Junior, iTanaka et all 
120111) . iFangetall d201ll) modeled the broadband emission of 
RCW 86 based on a semi-analytical solution to the NLDSA the- 
ory at SNR shocks. Their predictions are in conflict with the 
Fewra'-LAT upper limits. It should be noted that the SC spec- 
trum from electrons with the same index F = 1.8 is inconsistent 
with the radio spectral index deduced from the two flux points in 
Figure [4] To reproduce the radio data from MOST and Parkes, 
the magnetic field should be high enough to produce signifi- 
cant energy losses during the lifetime of the remnant (~ 1820 
y) and create a break in the SC spectrum. A reasonable fit is ob- 
tained by using a total energy injected into accelerated hadrons 
of ~ 7 x 10 49 d2 5 /n C m-3 erg and an average magnetic field of 
400 yt/G, with a break in the electron spectrum at ~ 40 GeV This 
magnetic field is much stronger than any estimates published so 
far on RCW 86, as discussed in section l4~2l and K ep ~ 2 x 10~ 3 
is much lower than that measured in the CR spectrum. 



4.3.2. Two-zone model 

The constraint on the magnetic field is clearly relaxed in a two- 
zone model in which the radio emission is produced by a sepa- 
rate electron population. In this case, a reasonable fit is shown 
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Fig. 4. Spectral energy distribution of RCW 86 with th e best-fit leptonic (to p) and hadronic (bottom) models. The radio data points 
are from Molonglo at 408 MHz and Parkes at 5 GHz dCaswell et alJ 1975L D. A. Green 2011, priv ate communication), and lower 
limits from MOST at 843 MHz and ATCA at 1.43 GHz dWhiteoak & Greenll996LlDickel et al.ll200ll) . X-ray spectra from ASCA and 
RXTE were obtained as des cribed in section[3l The .Fermz-LAT upper limits (99.9% C.L.) derived in section l2T2l and the H.E.S.S. data 
points in the VHE domain (lAharonian et al. 2009]) are also shown. For each plot, the blue (respectively red) lines denote the total 
broadband emission from the one-zone (respectively two-zone) modeling discussed in section|4] with SC, IC (with CMB only) and 
7T° spectra shown as solid, dashed, and dash-dotted lines, respectively. The associated best-fit parameters are provided in Table [2] 
Note that in the bottom panel, n a emission is similar for the two models, and equals the total gamma-ray emission in the one-zone 
model (solid blue line). In addition, the dotted black line shows the (excluded) case of n° emission from a E~ 2 hadron spectrum with 
T] p -0.15 d| 5 /n cm _ 3 and E maXiP = 100 TeV. 
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Table 2. Best-fit parameters of the leptonic and hadronic scenarios for the one-zone and two-zone models applied to the 
RCW 86 broadband emission (see Fig. ffli. N umbers in parentheses correspond to the constraints obtained when considering the 
Galactic ISRFs from iPorter & Strond d2005l) in the calculation of IC emission. Note that we do not consider any population of 
accelerated hadrons associated with the radio-SC-emitting electrons in the two-zone models (K ep is then left undetermined). We 
also fixed E maxp to 100 TeV in the leptonic scenarios. 



8 



Lemoine-Goumard et al.: The SNR RCW 86 using multi- wavelength observations 



in Figure [4] (bottom, red lines). It assumes a (lower limit on the) 
magnetic field of 50 //G in the X-ray-emitting zone (to suppress 
the IC component), with an injection spectrum with a power- 
law index of 1.8, an energy break at 3 TeV (produced by SC 
cooling), and an exponential cut-off at 20 TeV. This two-zone 
model does not modify the estimated energy injected into ac- 
celerated hadrons, and the spectrum of radio-emitting electrons 
is also unchanged with respect to the two-zone leptonic model 
( 14.2.2b . The resulting K ep is still much lower than 10~ 2 . In regard 
to the sp ectral shape deduc ed from these hadronic scenarios, we 
note that lVink et al.l d2006l) found that a concave lepton spectrum 
is favored over a simple power-law to explain the broadband SC 
emission from the NE region, even though they could not distin- 
guish a flattening at high energies T = 2.2 — » 2 from a flattening 
T = 2.2 — ¥ 1.5. The resulting 7r°-decay emission from such a con- 
cave hadron spectrum (given that relativistic lepton and hadron 
indices should be very similar) would certainly mimic an E~' 8 
spectrum at the highest energies, although it is unclear whether 
the lower-energy part of such a spectrum would still agree with 
the fermz-LAT upper limits. It is also interesting to note that, in 
a hadronic scenario, the gamma-ray emission is expected to be 
approximately proportional to the local density. Since the NW 
and SW regions both exhibit higher densities than the NE re- 
gion, these regions should be brighter in the TeV regime, assum- 
ing a uniform CR density along the SNR shell. However, while 
the H.E.S.S. image looks similar to th at observed in radio and 
X-rays, the azimuthal profile shown in lAharonian et al. 1 (120091) 
is well-fit with a constant, without any clear variation along the 
shell. Deeper VHE observations are required to investigate in 
more detail the gamma-ray morphology of RCW 86 and help 
discriminate the origin of the gamma-ray emission. 



5. Conclusion 

The Fermi-hPC[ upper limits in the HE domain derived in this 
work for the young SNR RCW 86 provide strong constraints on 
the injection spectrum of the primary population responsible for 
the extended VHE emission. A hadronic scenario can only re- 
produce the multi-wavelength data using a hard proton spectrum 
(spectral index T < 1.8) and a total energy injected into hadrons 
residing in RCW 86 of ~ 7 x 10 49 (>i/l cirT 3 )- 1 d z 25 erg. Given 
that the one-zone hadronic model suffers from several limitations 
(incompatible radio spectral index, very low K ep and extremely 
high B-field), only a two-zone model can fulfill all the observa- 
tional constraints, though still with a low K ep . In the one- and 
two-zone leptonic scenarios, the multi-wavelength data can be 
closely reproduced using electron spectral indices of ~ 2.0-2.3, 
a total energy injected in electrons of ~ 2 x 10 49 erg, and a rea- 
sonable average magnetic field of 15-25 /yG. In such a case, the 
most conservative upper limit to the total energy injected into 



2.5 



hadrons inside RCW 86 amounts to ~ 4 x 10 4y («/l cm 3 ) 1 d^ 
erg. 

These estimates of the total CR energy content in 
RCW 86 can be translated into CR pressure P t ot,CR = Et^ cR^y, 
where V is the effective CR volume within the SNR. Assuming 
that this volume is given by the best-fit parameters of the shell 
morphology observed with H.E.S.S. ( though not statistically 
significant over a uniform sphere, see lAharonian et al.l 120091) . 
we obtain V — (6 ± 2) x 10 59 d^ 5 cm 3 , and P t ot,CR must be 
conservatively lower than 3.8^° x 10~ 10 (E tot , p /7 x 10 49 erg) 
(n/0.1 cm" 3 )" 1 d" 1 , erg cm . This CR pressure, derived from 
the modeling of the broadband nonthermal emission of the en- 
tire SNR (see section 0), is very close to the CR pressure found 



bvlHelderet al.l (120091) in the NE region, and later confirmed by 
IVinketal.l(l2010h . Pmepp > 3.7 x 10~ 10 (n/0.1 cm 3 ) (kT/2.3 
keV) erg cirT 3 . If we were to apply the same comparison for 
the other regions along the RCW 86 shell, which all exhibit 
higher plasma densiti es (0.5-2 cm -3 , see lYamaguchi et ai1l201 U 
IWilliams et ail l20lTb . our up per limit to - Ptot.CR would most 
likely be inconsistent with the iHelder et al.l (120091) estimate re- 
garding the fractional CR pressure. However, the large uncer- 
tainties in the effective CR volume in the above calculations, 
among others, prevent us from drawing firm conclusions about 
the acceleration efficiency in RCW 86. 
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